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Abstract 

We report the discovery of a "dark area theorem," a new quantum optical 
relation for propagation of unmatched pulses in thick three-level A-type me- 
dia. We define dark area and derive the dark area theorem for a coherently 
prepared and inhomogeneously broadened lambda medium. We also obtain 
the first equation for the spatial evolution of the dark state amplitude prior 
to pulse-matching. 

Many advances have been made recently in controlling quantum systems with a pair 
of near-resonant optical pulses under conditions of relatively long-preserved two-photon co- 
herence among three quantum states. Key parameters being controlled include laser pulse 
intensity, frequency, duration, shape, timing and detuning. Well-known techniques and phe- 
nomena such as Rabi splitting |T|, spectroscopic dark-states 0, pulse matching ||, and 
anti-intuitive excitation have been exploited, with consequences including EIT (electro- 
magnetically induced transparency f5|), LWI (lasing without inversion ||), state-selective 
molecular excitation (STIRAP 0), and demonstrations of slow light || and fast light ||, 
to mention only a few prominent examples. 

However, in practically every instance a strong field is used to control a near-resonant 
medium so that a weak field can penetrate or amplify or excite transitions in a way normally 
forbidden. These are thin absorber applications, in the sense that the strong field's evolution 
in the medium is ignored. 

The understanding of pulse evolution in two-level media is much more complete. The 



McCall-Hahn Area Theorem |T0[ is available to provide a unifying picture of the way the 
medium acts back self-consistently. Even in situations where coherence is not complete, the 
two-level Area Theorem is a guide to the key role of self-consistent back action, and may 



suggest inverse applications in which the medium controls the field |TT|], i.e. thick absorber 
applications. 

In the present note we report the discovery of a new nonlinear propagation law for 
unmatched pulses in a three-level medium. This new propagation law, which is appropriately 
labelled a Dark Area Theorem, takes the form: 
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where the constant <ft an d the "dark area" 9 B are defined below. Here a and a 1 are the 
conventional Beer's Law absorption coefficients for the two transitions shown in Fig. ffl(a). 

The fundamental equations, from which the dark area theorem is derived, are familiar. 
For the physical pulses we have: 
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where the notation is indicated in Fig. |l|(a). The f2's are the two Rabi frequencies of 
the (assumed unchirped) fields: Q g = 2d g £ g /h, etc. The angular brackets are defined by 
((. . .)) = J F(A)(. . .)dA, which is an average over detunings arising from inhomogeneous 
broadening, where F is a normalized distribution, taken symmetric and Gaussian in the 
examples calculated below, with Fq setting the time scale for inhomogeneous relaxation. We 
take the same F for the two transitions. The \x parameters are related to the a's in the 
standard way for inhomogeneously broadened media: a = nfiF Q , and a' = 7r/i'F , and we 
do not assume \x = //. Finally, we have defined local-time coordinates r and ( in the frame 
propagating with velocity c in the medium: cr = ct — z and £ = z. 

In addition to the two field equations there are three equations for the atomic amplitudes 
g, g' and a. We prefer to write these equations in the dressed "bright-dark" basis sketched 
in Fig. |(b): 
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where we have defined the bright Rabi frequency as Q B = + anc i the dark Rabi 

frequency via the Fleischhauer-Manka relation |Q: Q B = —i(Q g Q g t — Q g tQ g )/Q 2 B . In these 
equations and throughout the paper the symbols B and D denote the bright and dark state 
amplitudes, which are related by a well-known field-dependent dressing rotation [13] in g - 
g' space: 
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where cos 9/2 = Q g /Q B and sin#/2 = Q g >/Q B - 

However, we have discovered that the ±2n7t ambiguity of this definition of 9 has physical 
meaning, and that the meaning is related to the spatial evolution of 9. This is a double 
surprise - that the ambiguity is physical and that the dressing angle should be treated as a 
propagation variable, actually as an "area" corresponding to the dark Rabi frequency: 
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where tan 2@((i ~ 00 ) = ^V(C> — °°)/^s(Cj —00). Now we explain these remarks by examin- 
ing the previously unexplored spatial evolution of 8. 
The first Dark Area propagation equation is: 

an nj 

^ = ^ + (aD*) - ^(aB*)) , (9) 

where the coefficients are ^-dependent and defined as fi + = + fi') + — fi) cos 8 and 
yU_ = — //) sin 8. This equation is derived from the definition of dark area and equations 
(§) and (fD, and with it we obtain the key to the propagation regime. An expression for 
(aD*) and (aB*) can be obtained in two steps, first by formal integration of the a equation 
up to a time T following the passage of the pulses, to obtain 

a(T; A) = j C e'^^ ({l B B)dr } (10) 

and second by substituting this integral into the right side of (|9[), where the resulting dou- 



ble integrals (over r and A) can be evaluated This is done in the regime of rapid 

inhomogeneous relaxation, in which Fq is a shorter time than the pulse durations or the 
medium's homogeneous lifetimes, and we obtain at time T: i(aD*) w — ^F q QbDoB and 
i{aB*) w — t;F QbBq. The subscript here denotes A = 0. 

Now we restrict attention to the case /1 = fi' for simplicity in describing the developments, 
in which case the central propagation equation (||) reduces to: 

98 D 2ifi 
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where 8b is the total dark area, i.e., 8d = 8((,t > T). Within the r > T regime we 
apply the solutions for B and D that are compatible with the asymptotic preparation of the 
medium, namely B = sin \ {8 D + 0) and D = cos^^ + 0), where tan(0/2) = g/g f taken 
at r — > —00. This converts (|TT|) directly into the Dark Area Theorem: 

-^ = --sm(0 D + 0), (12) 

which is easily seen to be the reduced form of ([]]) when a = a'. The more general equation 
(H) can be derived in the same way, simply by retaining a 7^ a' throughout. The Dark Area 
theorem makes asymptotic predictions for the pulses, and assigns distinct physical content 
to each of the branches of its solution shown in Fig. |[ 

As in the two-level case, the three-level area theorem does not predict pulse shapes, and 
an infinite variety of shapes can be asymptotically stable. In Fig. |3] we show the evolution 
of two initially unmatched and temporally offset pulses, with initial areas 1.57T and 4.87T. 
Propagation over 20 Beer's lengths shows highly non-trivial shifting and reshaping. In 
addition one can check that the final matched pulses conform to the dark area theorem's 
asymptotic rule: their amplitude ratio is given by tan(0/2). 
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A further surprise arising from the use of the dressing angle as a propagation variable 
is that the spatial behavior of the bright and dark amplitudes is already specified by the 
dressing transformation itself. In the r > T > F regime, where the space-time behavior of B 
and D is a self-consistent reponse to the pulses, we can write dB/d( = (dB / d9 D ) (d9n / d() , 
and a similar equation for D . With the aid of (|TT|), we can thereby obtain a new equation 
for D 2 alone: dD 2 /d( = aD 2 (l — D 2 ), with the elementary solution: 

D 2 (l - D 2 ) 

to (T±« = »«-«■ (13) 

where Di n here means the value of D at the arbitrary position of incidence Q n . The new D 2 
equation itself already shows that the alternate solution D = is not stable, which is exactly 
the conclusion reached in our earlier numerical study of the distinction between S7T-type 
and EIT-type propagation JT5] . 



More important, the analytic solution given here also provides the first compact pre- 
dictive expression |TB[ for spatial evolution of the dark state toward matched pulses, and 



confirms what has been known empirically, that the real pulses finally must become matched, 
which is the same as the limit D 2 — > 1. The straight lines of data points given in the inset 
in Fig. [| show the agreement between a completely numerical solution obtained from the 
original equation set @ - (H) and the simple expression (|13"1) predicted by dark area theory. 
The two different lines of points emphasize that for different times t/T* in the pulse the 
spatial evolution of \D\ 2 will attain the predicted asymptotic straight-line form more or less 
quickly. 

To summarize, we have presented the explicit form of an area theorem that is, to the 
best of our knowledge, the first to be found for propagation in thick three-level media. The 
most important ingredient of the derivation was the discovery that 9d, twice the dark-state 
dressing angle, is the key propagation variable for physical pulses in coherently prepared 



three-level media [|17]]. We have demonstrated that, in the domain of rapid inhomogeneous 
relaxation, approximate expressions for the response of the medium can be found, and that 
they lead directly to compact expressions for the spatial evolution of bright and dark state 
amplitudes, as in eqn. ([13|) , which were previously unknown. Numerical solutions show that 
these simple expressions are highly accurate. 

We expect applications of the dark area theorem to be numerous and interesting. One 
can show that the singular points of both tanOo and cot Op determine the number of zero- 
crossings of the physical pulse envelopes Q g and Q g ', such as the two that appear spon- 
taneously in Fig. § around aL = 5 and 10. Such zero-crossings appear closely related to 



phase jumps previously observed in Raman solitons [0. The inclusion of inhomogeneous 
broadening makes the dark area theorem capable of addressing three-level echo effects |L9 . 
Application of the singularity rules mentioned above will lead to non-trivial analogs of the 
McCall-Hahn rule for split-up of 2nn pulses. In conclusion, we mention again that the full 
expression of the theorem in eqn. allows fi ^ fjf. A detailed examination of ([I]) must 



be given elsewhere |[20|| , along with the spatial development of the bright field Qb and other 
elements missing here for lack of space. 
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FIGURE CAPTIONS 



Fig. [I] 

a) Excitation of the bare-state atom by two pulsed fields, denoted by their Rabi frequencies 
Q g and Q g i, and b) excitation of the counterpart dressed-state atom by bright and dark 
fields, as defined in the text below eqn. @. 

Fig. | 

Evolution of the dark area 9d with distance, as predicted by eqn. ([!]). Asymptotes are 
determined by 9d + 4> = 2n7r and tan(#o/2) tan(</>/2) = a/a'. The values used here are 
= 2t[/3 and a' = 4a. 

Fig. | 

Spatial snapshots of the evolution of two initially unmatched pulses with areas 1.57T and 
4.87T, through 20 Beer's lengths of a A medium with = 27r/3, obtained by exact numerical 
solution of eqns. as a test of eqn. fll2|) . The inhomogeneous detuning distribution 

F(A) is taken Gaussian, and local time is shown in units of T* = \^2nF . At aC, = 20 the 
pulses are almost matched, with Q g i/Q g « 1.6, in good agreement with the final amplitude 
ratio predicted by the branch value of 9n, i.e., tan(27r — 4>) = tan(4-7r/3) = \^3. 

Fig. 1 

Dark state population as a function of time is shown for five different depths of propagation 
of the same pulses shown in Fig. |3|. Evolution to \D\ 2 = 1 is evident. The specific prediction 
of eqn. (0) is also checked by computing \D((,t)\ 2 for these pulses over 8 equal intervals 
at a( = 0, 2.5, 5.0, . . . , 20.0 at both times r/T* = 15 and 30, indicated by vertical dashed 
lines in the right-hand half. The two sets of straight-line data with approximately unit slope 
shown in the inset provide excellent confirmation of the dark area theory. 
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FIG. 2. J.H. Eberly and V.V. Kozlov, "Dark Area Theorem" 
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FIG. 3. J.H. Eberly and V.V. Kozlov, "Dark Area Theorem" 
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FIG. 4. J.H. Eberly and V.V. Kozlov, "Dark Area Theorem' 
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